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We report progress on a non-perturbative calculation of the light quark mass renormalization 
factor Zm S,LAT (n,a), using dynamical Asqtad fermions. This quantity is used to determine the 
light quark masses in the conventional MS scheme. Such a non-perturbative determination of 
Z,„ removes uncertainties due to truncation of its perturbative expansion, currently known to two 
loops. These calculations have been carried out using publicly available MILC lattices with 
spacings of w 0.125 and 0.09 fm. 
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1. Introduction 

Quark masses are fundamental parameters of the standard model, and as such are inputs whose 
determination must ultimately come from experimental measurement. Lattice simulations of QCD 
provide a means by which these parameters may be precisely determined and their values systemat- 
ically improved upon. Because quark masses are scheme and scale dependent, one must determine 
"Z-factors" to convert the bare lattice masses determined from simulations to a continuum scheme. 

Z-factors may be determined using lattice perturbation theory, but these determinations nec- 
essarily suffer from uncertainties due to truncation of the perturbative expansion. Extending these 
calculations to higher order is difficult, particularly when using improved actions as is done in most 
current simulations. Another strategy, more amenable to systematic improvement, is to compute 
Z-factors non-perturbatively using the lattice simulation itself. 

In this article we report on progress in this area using staggered quarks and the Asqtad [1] [2] [3] 
O (a 4 , a 2 Ok) -improved action. This is the first such unquenched study using staggered quarks. First 
we describe the methodology of non-perturbative renormalization (NPR) [4] used in the study, 
highlighting details specific to the use of staggered fermions. Then we present our results, and 
conclude with a discussion of further work to be done. 



2. Methodology 
2.1 Basic principles 

Quark masses are scheme and scale dependent quantities. For a given lattice regularization, 
one tunes the bare input masses to the physical point of QCD. To compare amongst different deter- 
minations, these quantities must be converted to a common scale and scheme. This is convention- 
ally taken to be the MS scheme at a scale of jU = 2GeV: 

m m ^)=Z^{^a) a -^. (2.1) 

The factor Zm S,LAT (jj.,a) has been calculated for the Asqtad action to two loops using lattice per- 
turbation theory [5]. 

Because the MS scheme is a continuum scheme, one must first convert to a scheme that is sen- 
sible in both lattice and continuum theories, and then use an additional matching factor to convert 
between schemes. We take the intermediate scheme to be the RI/MOM scheme l , to be explained 
momentarily: 

zf' Mr (/M) = ZF^QI) U R \^p) Z^(p,a). (2.2) 

Here U RI {jJL,p) is a perturbative running factor computed in the continuum, known to 4 loops, and 
Zm S S/ (aO is the MS to RI conversion factor, calculated to 3 loops in the continuum [7]. If one was 
content to work in the RI scheme, this latter factor could be avoided entirely. 

In the RI/MOM scheme, Z-factors are defined such that Green's functions in a fixed (Landau) 
gauge, suitably projected, take their tree-level values. Such Z-factors may be computed directly 



'in practice we use the RI' scheme. This differs from the RI scheme in the definition of the wavefunction renormal- 
ization factor X q (2.3) and results in a difference of 0(1%) [6]. 
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from lattice simulations [4]. One only requires that the scale of external momenta p is chosen such 
that Aqcd <C \p\ <C \ r This ensures that non-perturbative contributions to Z-factors and lattice 
artifacts are simultaneously under control. Thus the NPR method trades the truncation uncertainty 
of a perturbative evaluation for uncertainties due to lattice artifacts, non-perturbative contributions, 
and statistics. These, however, are more amenable to systematic reduction. 

In principle, one could extract Z,„ using only the following projections of the quark propagator: 



Z q (p) 



1 



1 



-Tr 



=Z q (p) 



(y^oi^tp] 



Z m [p)m + C\— T - + ... 



(2.3) 



(2.4) 



The propagator S\p) is a 16-by-16 matrix in spin-taste space, as are the double-barred projection 
matrices. The latter are the staggered analogues of their continuum counterparts. In practice, 
the projection (2.4) suffers severely from non-perturbative contamination due to the condensate 
term [8]. One must make recourse to Ward identities to obviate this issue. 

The lattice Ward identities of interest relate the scalar and pseudoscalar bilinear operators to 
the quark propagator: 
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(2.5) 



(2.6) 



Here the A are amputated vertices of the bilinear operator with the indicated spin-taste structure. 



In perturbation theory, one has the identity Z„ 



_j_ 

z s 



-. On the lattice, the equality is spoiled 



by non-perturbative contributions, which can be seen from inserting (2.4) into (2.5) and (2.6). One 
expects that the scalar bilinear in (2.6) gives a closer approximation to Z m because the derivative 



dm 



will eliminate the \ condensate pole. 



2.2 Implementation 

Computations are being earned out using the Chroma [9] software library for Lattice QCD. 
We have added support for staggered momentum-source inversions and written routines to compute 
the bilinears of interest. All configurations are first fixed to Landau gauge. For a given momentum 
p' in the reduced Brillouin zone (0 < pL < 16 inversions are carried out on momentum sources 
to produce the correlators: 



<X(x)X(-p' + XB) > = <x{x)$ B (-p') > . 



(2.7) 



Here B is one of the sixteen binary vectors e.g. (1, 1,0,0). Note that these momentum sources, 
being simple phase factors of the form e' px , differ from those used in the quenched study [10], 
and result in a propagator that does not explicitly break taste symmetry. From these correlators 
the "polespace" propagators are constructed by performing sixteen Fourier transforms on the free 
space index: 

Sab\p'] = [ d 4 xeW + " A > x < %{x)U-p) > • (2-8) 
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Thus the propagators are 16-by-16 matrices in spin- taste space. 

Inversions are carried out on publicly available MILC lattices [11] with 2+1 flavors of Asqtad 
quarks, and are computed on Fermilab clusters using resources from a USQCD grant. We use both 
fine (a m 0.09 fm) and coarse (a « 0.125 fm) lattices. On the coarse lattices, we have varied both 
the sea and valence masses to produce unquenched data with masses am = .03, .02, .01. On the 
fine lattices we have studied masses am va i = .018, .012, .006 and am sea = .0062, and intend to alter 
the sea quark masses to produce unquenched data. 

We study the propagators and bilinears over a range of momenta to understand the non- 
perturbative and discretization effects. The momenta are chosen using the naive "cylindrical cut" 
method [12] to reduce discretization effects. Data shown in this paper are all obtained using 8 
lattices. This results in relatively small statistical uncertainties because of our use of momentum 
sources. 

3. Results 

Figure 1 shows typical results for the scalar and pseudoscalar bilinears. The l/p 2 non- 
perturbative contamination is clear in the graph of Zp 1 , but is absent in Z^T 1 , consistent with the dis- 
cussion in section 2. Figure 2 shows a closeup of the same Z^T 1 data, as well as the "scale-invariant" 
result, ZfJ, obtained using the perturbative running factor to run the data to a common scale, here 
(ap) 2 = 2. If discretization effects were negligible, one would expect to be a constant, at least 
in the region of large (ap) 2 where the perturbative expansion is sensible and non-perturbative con- 
tributions are small. Discretization effects proportional to (ap) 2 would result in the curve being a 
straight line with non-zero slope. The observed data is consistent with this. Thus we remove arti- 
facts by using the intercept of the interpolated line. This is detailed in Figure 3. The fitted points 
are chosen to minimize the % 2 per dof of a straight line fit. In Figure 4 we have plotted the result of 
such a procedure for the three masses on the coarse MILC lattices. Also shown is a naive linear fit 
to the data, which, when combined with the bare strange mass and lattice spacing as determined by 
the MILC collaboration [11], gives a final strange quark mass of mf s = \06(6)MeV on the coarse 
lattices. This is to be compared with the perturbative value of mf s = S4(5)MeV obtained on the 
coarse lattices by the MILC collaboration [5]. 

4. Conclusions, Outlook, Future Work 

It should be noted that the above analysis is preliminary. A linear chiral extrapolation has been 
carried out, but a glance at Figure 4 shows that statistics must be improved to verify whether this is 
permissible. In addition, analogous calculations must be performed on the fine lattices to permit a 
continuum extrapolation. We also intend to calculate Z q using the vector and axial-vector bilinears, 
as this makes use of the RI scheme directly and has no uncertainty associated with the choice of 
definition of Z q [13]. 

Current data on the coarse lattices indicate a strange quark mass significantly larger than that 
obtained by perturbative analysis. The fine lattice data also appear to support this result, though a 
fully unquenched analysis is pending. This seems to be consistent with a general trend that the non- 
perturbatively obtained masses are larger than those obtained perturbatively [14] [15]. This could 
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Figure 1: Z p and Z s on the MILC coarse lattice ensemble (8 configs). The \/p non-perturbative con- 
tamination in 1 is clearly visible. 



N 



L.S 



a~.l25fm 
am = .0 1 



rySI 



(ap) 2 

Figure 2: Result for Z m = Z^T 1 . After applying the perturbative running factor to obtain Z%, the data 
noticeably flattens out and is consistent with a straight line at higher (ap) 2 . 
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be due to deficiencies of the perturbative expansions, or to systematic non-perturbative effects that 
have not been accounted for. The use of non-exceptional momenta should help to distinguish these 
possibilities [15]. 
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